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Abstract 

The extragalactic y-ray sky is dominated by the emission arising from blazars, one of the most pecu¬ 
liar classes of radio-loud active galaxies. Since the launch of Fermi several methods were developed to 
search for blazars as potential counterparts of unidentified y-ray sources (UGSs). To confirm the nature 
of the selected candidates, optical spectroscopic observations are necessary. In 2013 we started a spec¬ 
troscopic campaign to investigate y-ray blazar candidates selected according to different procedures. The 
main goals of our campaign are: 1) to confirm the nature of these candidates, and 2) whenever possible 
determine their redshifts. Optical spectroscopic observations will also permit us to verify the robustness 
of the proposed associations and check for the presence of possible source class contaminants to our 
counterpart selection. This paper reports the results of observations carried out in 2014 in the Northern 
hemisphere with Kitt Peak National Observatory (KPNO) and in the Southern hemisphere with the South¬ 
ern Astrophysical Research (SOAR) telescopes. We also report three sources observed with the Magellan 
and Palomar telescopes. Our selection of blazar-like sources that could be potential counterparts of UGSs 
is based on their peculiar IR colors and on their combination with radio observations both at high and 
low frequencies (i.e., above and below ~ 1 GHz) in publicly available large radio surveys. We present the 
optical spectra of 27 objects. We confirm the blazar-like nature of 9 sources that appear to be potential 
low-energy counterparts of UGSs. Then we present new spectroscopic observations of 10 active galaxies 
of uncertain type associated with Fermi sources, classifying all of them as blazars. In addition, we present 
the spectra for five known y-ray blazars with uncertain redshift estimates and three BL Lac candidates 
that were observed during our campaign. We also report the case for WISE J173052.85-035247.2, can¬ 
didate counterpart of the source 2FGL J1730.6-0353, which has no radio counterpart in the major radio 
surveys. We confirm that our selection of y-ray blazars candidates can successfully indentify low-energy 
counterparts to Fermi unassociated sources and allow us to discover new blazars. 

Subject headings: galaxies: active - galaxies: BL Lacertae objects - radiation mechanisms: non-thermal 

1. Introduction 

Blazars are the largest population of active galactic 
nuclei (AGNs) detected in the y range (Abdo et al. 2010; 
Nolan et al. 2012). Their non-thermal emission ex- 
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tends from radio to TeV energies and it is coupled 
with very rapid variability, high and variable polar¬ 
ization, superluminal motion, high luminosities (e.g., 
Urry & Padovani 1995) and peculiar infrared (IR) col¬ 
ors (Massaro et al. 2011a). Since 1978, well before 
the establishment of the unification scenario for the 
AGNs, their peculiar properties were described in 
terms of emission arising from particles accelerated 
in a relativistic jet closely pointed along our line of 
sight (Blandford & Rees 1978). 

Adopting the nomenclature of the blazar subclasses 
described in the Roma-BZCAT: Multi-frequency Cata- 


these selection procedures of y-ray blazar candidates 
(e.g.. Stem & Assef 2013) and such degeneracy can 
be only removed with optical spectroscopic obser¬ 
vations (e.g., Masetti et al. 2013; Shaw et al. 2013a; 

Shaw et al. 2013b; Paggi et al. 2014; Massaro et al. 2014a). 
Furthermore, a detailed knowledge of the num¬ 
ber of UGSs is of paramount importance for in¬ 
stance to provide constraints on dark matter models 
(Abdo et al. 2014). Many UGSs could be blazars, but 
how many of them are actually blazars is still unknown 
due in part to the incompleteness of catalogs used for 
the associations (Ackermann et al. 2011a). 


logue of Blazars^ (e.g., Massaro et al. 2009; Massaro et al. 2011b),Thus, motivated by these arguments, we started an 


we distinguish between BL Lac objects (i.e., BZBs) 
and the fiat spectrum radio quasars indicated as BZQs. 

The former present optical spectra with emission 
and/or absorption lines of rest frame equivalent width 
EW < 5 A (e.g., Stickel et al. 1991; Stoke et al. 1991; 
Laurent-Muehleisen et al. 1999) while the latter show 
typical quasar-like optical spectra with strong and 
broad emission lines. In the Roma-BZCAT there are 
also several BZBs indicated as BL Lac candidates; 
these sources were indicated as BL Lacs in literature 
and thus reported in the catalog but, lacking their op¬ 
tical spectra, a correct classification is still uncertain 
(see also Massaro et al. 2014a). 

We developed several methods to search for y-ray 
blazar candidates that could be counterparts of the 
unidentified y-ray sources (UGSs, 2FGL Nolan et al. 2012) 
on the basis of the peculiar IR colors of known y-ray 
blazars (e.g. D’Abrusco et al. 2013; Massaro et al. 2013a; 
Massaro et al. 2013b), discovered thanks to all-sky 
survey of the Wide-Field Infrared Survey Explorer 
(WISE; Wright et al. 2010), or using radio observa¬ 
tions in combination with IR colors (Massaro et al. 2012c) 
as well as at low frequencies in the MHz regime 
(Massaro et al. 2012b; Nori et al. 2014). In addition, 
multifrequency analysis based on X-ray follow-up ob¬ 
servations (e.g., Mirabal & Halpern 2009; Paggi et al. 2013; 
Takeuchi et al. 2013; Stroh et al. 2013; Acero et al. 2013) 
and radio campaigns (e.g., Petrov et al. 2013) were 
performed to search for potential UGS counterparts. 

Since the positional uncertainty of Eermi is a few 
tenths of a degree, all the proposed methods and 
the multifrequency follow-up observations are pri¬ 
marily useful to decrease the number of potential 
counterparts for the UGSs. There could be pos¬ 
sible contamination by different source classes in 


http://www.asdc.asi.it/bzcat/ 


optical spectroscopic campaign in 2013 aiming to con¬ 
firm the real nature of the proposed low-energy coun¬ 
terparts of UGSs selected according to our methods. 

In this paper we report the results of observations 
carried out since January 2014 in both hemispheres. 
We mainly used the Kitt Peak National Observatory 
(KPNO) for our targets in the Northern hemisphere in 
addition to two more observations performed at Palo- 
mar. For targets mainly visible in the Southern sky, 
we present spectra obtained with the Southern Astro- 
physical Research (SOAR) telescope and one more ob¬ 
servation carried out at Magellan. Preliminary results 
for our exploratory program in the Northern hemi¬ 
sphere obtained with the Telescopio Nazionale Galileo 
(TNG), the Multiple Mirror Telescope (MMT) and the 
Observatorio Astronomico Nacional (OAN) in San Pe¬ 
dro Martir (Mexico) were already presented in Paggi 
et al. (2014). In addition, the results of our obser¬ 
vations carried out in the 2013 campaign with SOAR 
and KPNO are described by Landoni et al. (2014) and 
Massaro et al. (2014b), respectively. 

The paper is organised as follows: Section 2 con¬ 
tains the different methods that we employed for the 
sample selection, in Section 3 we present our observa¬ 
tional setup and we discuss the data reduction proce¬ 
dures. Then in Section 4 we describe the results of our 
analysis for different types of sources in our sample. 
Finally, Section 5 is devoted to summary and conclu¬ 
sions. We use cgs units unless otherwise stated. 

2. Selecting the sample 

The surveys and the catalogs used to search for the 
counterparts of our targets are listed in the following. 
These symbols are also reported in Table 1. Below 
IGHz we used the VLA Low-Frequency Sky Survey 
Discrete Source Catalog (VLSS; Cohen et al. 2007, - 
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V) and the recent revision VLLSr^ (Lane et al. 2014), 
the Westerbork Northern Sky Survey (WENSS; Rengelink et al. 

- W), the Sydney University Molonglo Sky Survey 
(SUMSS; Mauch et al. 2003, - S), the Parkes-MIT- 
NRAO Surveys (PMN; Wright etal. 1994, - Pm), 
the Parkes Southern Radio Source catalog (PKS; 

Wright et al. 1990, - Pk), and the Low-frequency 
Radio Catalog of Flat-spectrum Sources (LORCAT; 

Massaro et al. 2014b, - L). At higher radio frequen¬ 
cies we also verified the counterparts in the NRAO 
VLA Sky Survey ^ (NVSS; Condon et al. 1998, - 
N), the Australia Telescope 20 GHz Survey (AT20G; 

Murphy et al. 2010, - A), the Combined Radio All- 

Sky Targeted Eight-GHz Survey (CRATES; Healey et al. 2007, 

- c). In the infrared, we queried the WISE all-sky sur¬ 
vey in the AllWISE Source catalog"^ (Wright et al. 2010, 

- w) and the Two Micron All Sky Survey (2MASS; 
Skrutskie et al. 2006, - M) since each WISE source 
is automatically matched to the closest 2MASS po¬ 
tential counterpart (see Cutri et al. 2012, for details). 

Then, we also searched for optical counterparts, with 
or without possible spectra available, in the Sloan 
Digital Sky Survey Data Release 9 (SDSS DR9; 
e.g. Ahn et al. 2012, - s), in the Six-degree-Field 
Galaxy Redshift Survey (6dFGS; Jones et al. 2004; 

Jones et al. 2009, - 6). At high-energies, in the 
X-rays, we searched the ROSAT all-sky survey in 
both the ROSAT Bright Source Catalog (RBSC; 

Voges et al. 1999, - X) and the ROSAT Faint Source 
Catalog (RFSC; Voges et al. 2000, - X), as well as 
XMM-Newton Slew Survey (XMMSL; Saxton et al. 2008; 
Warwick et al. 2012, - x), the Deep Swift X-Ray Tele¬ 
scope Point Source Catalog (IXSPS; Evans et al. 2014, 

- x), the Chandra Source Catalog (CSC; Evans et al. 2010, 

- x) and the Swift X-ray survey for all Eermi UGSs 
(Stroh et al. 2013). Note that we use the same symbol 
for the X-ray catalogs of XMM-Newton, Chandra and 
Swift. 

Our sample lists a total of 27 sources divided as 
follows: 9 are UGSs for which the analysis based on 
the IR colors of blazar candidates indicated them as 
blazar-like sources that were observed during our cam¬ 
paign; 10 are indeed classified as active galaxies of un¬ 
certain type (AGUs) according to The Second Catalog 
of Active Galactic Nuclei Detected by the Fermi Large 


^http://heasarc.gsfc.nasa.gov/W3Browse/all/vlssr.html 

^http://heasarc.gsfc.nasa.gov/W3Browse/radio-catalog/nvss.html 

^http://wise2.ipac.caltech.edu/docs/release/allwise/ 


Area Telescope (2LAC; Ackermann et al. 2011a) for 
no optical spectra were available when the cat¬ 
alog was released and have been observed as part 
of our campaign. Some sources have been also 
selected based on low-frequency radio information 
(Massaro et al. 2013b). The remaining 8 sources are 
BL Lac candidates and BL Lacs, either detected or 
not by Eermi listed in the Roma-BZCAT for which no 
optical spectroscopic information were found in litera¬ 
ture (Massaro et al. 201 lb) or with uncertain/unknown 
redshift estimate. 

We discuss our spectroscopic analysis in the follow¬ 
ing subsections while in Table 1 we summarize our 
results and report the multifrequency notes for each 
source with the only exceptions of those listed in the 
Roma-BZCAT. The finding charts for all the sources 
are shown in Figures 1-3 using the USNO-Bl Cata¬ 
log (Monet et al. 2003) and the Digitized Sky Survey^. 
We highlight that some of the sources observed during 
our campaign have also been observed at different ob¬ 
servatories and groups. However we re-observed these 
targets for two main reasons: 1) when our observations 
were scheduled and performed these spectra were not 
yet published; 2) due to the well-known BL Lac vari¬ 
ability in the optical energy range, there is always the 
chance to observe the source during a quiescent or low 
state and detect some emission and/or absorption fea¬ 
tures that could allow us to constrain their redshifts. 

3. Observations and Data Reduction 

Optical spectra of most of the sources accessi¬ 
ble from the Northern hemisphere were acquired in 
remote observing mode with the KPNO Mayall 4- 
m telescope using the R-C spectrograph, while the 
sources in the Southern hemisphere were observed 
remotely with the SOAR 4-m telescope using the 
High Throughput Goodman spectrograph. We also 
present the optical spectrum of the source WISE 
J024440.30-581954.5, observed on UT 2014 Febru¬ 
ary 3 in visitor mode at Las Campanas Observatory 
using the Magellan 6.5-m telescope in combination 
with the IMACS instrument. Finally, we present opti¬ 
cal spectra of WISE J025333.64-r321720.8 and WISE 
J085654.85-r714623.8 obtained in visitor mode on 
UT 2014 February 22 with the Double Spectrograph 
(DBSP) at the Hale 200-inch Telescope at Palomar Ob¬ 
servatory. Our scientific goal, i.e. the classification of 


^http://archive.eso.org/dss/dss 


3 




the selected targets, is best achieved with broad spec¬ 
tral coverage. We thus adopted a slit width of 1.2" (1") 
and the low resolution gratings yielding a dispersion of 
~3 (2) A pixel“^ for KPNO (SOAR). Our observations 
took place the nights UT 2014 June 4 and 5 at KPNO 
and on UT 2014 April 24 at SOAR during grey time. 
The average seeing for both runs was about 1" and 
conditions were clear. Additional observations were 
made on 2014 February 3 with the 6.5m Baade Magel¬ 
lan telescope using the Inamori Magellan Areal Cam¬ 
era and Spectrograph (IMACS; Bigelow et al. 1998). 
The f/2 camera was used in combination with the 3001 
mm“^ grism (blaze angle 17.5 degrees) and a 0.7" slit 
to yield spectra with dispersion of 1.34 A pixeU^ and 
FWHM resolution of ~4 A. Conditions were photo¬ 
metric and the seeing was generally < 1". The whole 
set of spectroscopic data acquired at these telescopes 
was optimally extracted (Horne 1986) and reduced fol¬ 
lowing standard procedures with IRAF (Tody 1986). 
For each acquisition we performed bias subtraction, 
flat fleld correction and cosmic rays rejection. Since 
for each target we secured two individual frames, we 
averaged them according to their signal to noise ratios 
(SNR). For questionable spectral features, we have ex¬ 
ploited the availability of two individual exposures in 
order to better reject spurious signals. The wavelength 
calibration was achieved using the spectra of He-Ne- 
Ar or Hg-Ar lamps that assure coverage of the entire 
range. To take into account drift and flexures of the 
instruments during the night, we took an arc frame 
before each target to guarantee a good wavelength 
solution for the scientiflc spectra. The achieved ac¬ 
curacy is about ~0.3 (0.4) A rms for KPNO (SOAR). 
Although our program did not require accurate pho¬ 
tometric precision, we observed a spectrophotometric 
standard star to perform relative flux calibration on 
each spectrum. Finally, we corrected the whole sam¬ 
ple for the Galactic absorption assuming Eb-v values 
computed using the Schlegel et al. (1998) and Cardelli 
et al. (1989) relations. To better detect faint spectral 
features and measure redshifts, we normalised each 
spectrum to its continuum. For the Palomar observa¬ 
tions, we observed the candidates through a 1.5" slit 
with the ~ 5500 A dichroic to split the light across 
the blue and red arms of DBSP. Both sides have had 
resolving powers R = T/A/l ~ 1000, and the data were 
reduced as above. 


4. Results 

4.1. Unidentified Gamma-ray Sources 

Here we provide details for the low-energy coun¬ 
terparts of the 9 UGSs observed in our sample. All 
these lT/5'£’-selected counterparts where found in 
the analysis performed by Massaro et al. (2013a) 
to have IR WISE colors similar to known y-ray 
blazars. All of them were predicted to be BZB-like 
sources having the IR colors more consistent with 
those of the Eermi BZBs rather than the BZQs (see 
D’Abrusco et al. 2013, for more details). The spec¬ 
tra of the whole UGSs listed in Table 1 are shown 
in Figure 4-12. The blazar counterpart for one of 
these sources, WISE J173052.85-035247.2 (candidate 
counterpart of 2FGL J1730.6-0353) appears intrigu¬ 
ing. On the basis of our optical spectra, we classify 
the source as a BL Lac object and we put a lower limit 
on its redshift of 0.776 based on the detection of an 
intervening doublet system of Mg II {AAobs = 4965 
- 4977 A with EWobs = 3.4 - 2.1 A; see also Fig¬ 
ure 9). However, this source is not associated with 
any NVSS sources, as expected for BL Lac objects. 
We note that, even if deeper radio observations de¬ 
tect emission from the source, blazars are traditionally 
defined as radio-loud sources on the basis of current 
radio data. All confirmed blazar in Roma-BZCAT are 
indeed detected at 1.4 GHz with fluxes above a few 
mJy, and radio-quiet blazars are extremely rare objects 
(Londish et al. 2004; Paggi et al. 2013). In particular, 
only 14 BZBs out of the 1220 present in the Roma- 
BZCAT have a radio flux density at L4GHz lower 
than 2.5 mJy that is the average flux limit of the NVSS 
survey^. 

There are two sources, WISE J174507.82-1-015442.4 
and WISE J174526.95-1-020532.6, potential counter¬ 
parts to 2FGL J1745.6-r0203 found on the basis of 
the IR color selection method. According to our 
multifrequency investigation we note that only WISE 
J174507.82-r015442.4 has a radio counterpart, thus we 
conclude that W75'£’ J174526.95-r020532.6 is a normal 
radio-quiet quasar contaminating our selection. How¬ 
ever we cannot firmly establish the real blazar nature of 
WISE J174507.82-1-015442.4 since the lack of multi¬ 
wavelength radio observations did not allow us to ver¬ 
ify the flatness of its radio spectrum. For both objects 
we determined their redshifts based on broad emission 
lines (Hof and H/5 are both present, see Figures 10-11). 


^http://heasarc.gsfc.nasa.gov/W3Browse/radio-catalog/nvss.html 
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Table 1: Selected sample and observation log. The sample is divided into three subsamples that are: potential coun¬ 
terparts of the UGSs, sources classified as AGUs, BL Lac candidates and BL Lacs with uncertain redshift listed 
in the Roma-BZCAT. Among the Roma-BZCAT sources, those marked with an asterisk are detected by Fermi 


(Nolan et al. 2012). Symbols used for the multifrequency notes are described in Sec. 2. 


Name 

Fermi 

name 

R.A. 

(J2000) 

Dec. 

(J2000) 

Obs. Date 
(yyyy-mm-dd) 

Exp. 

(s) 

notes 

class 

z 


J112325.37-252857.0 

2FGLJ1123.3-2527 

11:23:25.37 

-25:28:57.0 

2014-04-23 

300 

N,w,M,6,U,g - (z=0.145784-QSO-Jones+09) 

QSO 

0.148 

J125949.80-374858.1 

2FGLJ1259.8-3749 

12:59:49.80 

-37:48:58.2 

2014-04-23 

600 

S,N,w,U 

BL Lac 

? 

J132840.61-472749.2 

2FGLJ1328.5-4728 

13:28:40.43 

-47:27:48.7 

2014-04-23 

600 

S,rf,w,M,U,g,u,x - SED in Takeuchi+13 

BL Lac 

? 

J134042.01-041006.8 

2FGLJ 1340.5-0412 

13:40:42.01 

-04:10:07.2 

2014-06-05 

1200 

N,E,w,U,u,x - SED in Takeuchi+13 

BL Lac 

? 

J134706.88-295842.4 

2FGLJ 1347.0-2956 

13:47:06.88 

-29:58:42.5 

2014-04-23 

900 

S,N,rf,w,M,U 

BL Lac 

? 

J173052.85-035247.2 

2FGLJ1730.6-0353 

17:30:52.86 

-03:52:47.2 

2014-06-05 

900 

w,M 

BL Lac 

>0.776 

J174507.82+015442.4 

2FGLJ1745.6+0203 

17:45:07.82 

+01:54:42.6 

2014-06-05 

900 

N,w,M 

QSO 

0.078 

J174526.95+020532.6 

2FGLJ1745.6+0203 

17:45:26.96 

+02:05:32.8 

2014-06-05 

1200 

w,M 

QSO 

0.335 

J183535.34+134848.8 

2FGLJ1835.4+1349 

18:35:35.35 

+ 13:48:48.8 

2014-06-05 

600 

V,T,N,87,rf,w,M,U 

BL Lac 

7 


J025333.64+321720.8 

2FGLJ0253.4+3218 

02:53:33.64 

+32:17:20.9 

2014-02-22 

600 

N,87,GB,c,rf,w 

QSO 

0.859 

J074642.31-475455.0 

2FGL J0746.5-4758 

07:46:42.31 

-47:54:55.2 

2014-04-23 

600 

Pm,S,A,c,rf,w,M 

BL Lac 

7 

J084502.48-545808.4 

2FGL J0844.8-5459 

08:45:02.47 

-54:58:08.6 

2014-04-23 

600 

Pm,A,rf,U,X 

BL Lac 

7 

J085654.85+714623.8 

2FGLJ0856.0+7136 

08:56:54.86 

+71:46:23.9 

2014-02-22 

600 

W,N,87,GB,c,rf,w,M,g,X 

QSO 

0.542 

J114141.80-140754.6 

2FGLJ 1141.7-1404 

11:41:41.84 

-14:07:53.5 

2014-04-23 

900 

L,Pm,N,c,w,M,U 

BL Lac 

7 

J123824.39-195913.8 

2FGLJ1238.1-1953 

12:38:24.40 

-19:59:13.5 

2014-04-23 

900 

Pm,N,A,c,rf,w,M,g,X 

BL Lac 

7 

J140609.60-250809.2 

2FGLJ1406.2-2510 

14:06:09.60 

-25:08:09.3 

2014-04-23 

600 

L,N,w,M,U,u,x - SED in Takeuchi+13 

BL Lac 

7 

J162638.15-763855.5 

2FGLJ 1626.0-7636 

16:26:38.18 

-76:38:55.5 

2014-04-23 

300 

Pm,S,c,rf,w,M,g 

BL Lac 

0.1050 

J184931.74+274800.8 

2FGLJ1849.5+2744 

18:49:31.69 

+27:48:00.9 

2014-06-04 

1800 

W,N,87,c,rf,w,M 

BL Lac 

7 

J195945.66-472519.3 

2FGLJ1959.9-4727 

19:59:45.48 

-47:25:19.3 

2014-04-23 

600 

S,w,M,U,g,u,X,x - SED in Takeuchi+13 

BL Lac 

>0.519 


J024440.30-581954.5 

BZBJ0244-5819 

02:44:40.31 

-58:19:54.6 

2014-02-03 

600 

BL Lac candidate at z=0.265 

BL Lac 

>0.265 

J121752.08+300700.6 

BZBJ1217+3007* 

12:17:52.09 

+30:07:00.7 

2014-06-04 

300 

BL Lac at z=0.13? 

BL Lac 

7 

J 135949.71-374600.7 

BZBJl 359-3746* 

13:59:49.72 

-37:46:00.8 

2014-04-23 

300 

BL Lac 

BL Lac 

7 

J155333.56-311830.9 

BZBJ1553-3118* 

15:53:33.56 

-31:18:31.0 

2014-04-23 

300 

BL Lac (z=0.132-BLLac-Masetti+13) 

BL Lac 

7 

J164924.98+523515.0 

BZBJ1649+5235* 

16:49:25.00 

+52:35:15.0 

2014-06-05 

900 

BL Lac candidate at z=? 

BL Lac 

7 

J170209.63+264314.7 

BZBJ1702+2643 

17:02:09.64 

+26:43:14.8 

2014-06-05 

1200 

BL Lac candidate at z=? 

BL Lac 

7 

J180945.39+291019.8 

BZBJl 809+2910* 

18:09:45.39 

+29:10:20.0 

2014-06-05 

900 

BL Lac candidate at z=? 

BL Lac 

7 

J203923.51+521950.1 

BZBJ2039+5219* 

20:39:23.50 

+52:19:49.9 

2014-06-04 

600 

BL Lac candidate at z=0.053 

BL Lac 

7 


Their redshifts are reported in Table 1 . 

4.2. Gamma-ray Active Galaxies of Uncertain 
type 

In this subsection we discuss the AGUs in our ob¬ 
served sample. The multifrequency notes relative to 
each source are reported in Table 1, as previously 
done for the UGSs. The spectra of the whole AGU 
sample are shown in Figure 13-22. Our spectro¬ 
scopic observations confirm that 8 out of 10 sources 
are BL Lac objects while the remaining two, hav¬ 
ing quasar-like optical spectra and fiat radio spec¬ 
tra, are indeed BZQs. The optical spectra of these 
two BZQs, WISE J025333.64-r321720.8 and WISE 
J085654.85-r714623.8, contain broad emission lines 
identified as Mg II and HyS (see Figures 13 and 16 
for more details), and permit us to measure their red- 
shifts of 0.859 and 0.542, respectively. In particular, 
the identifications of [Oil], HyS and [OIII] doublet in 
the optical spectrum of WISE J025333.64-r321720.8 
are uncertain due to low SNR. The optical spectrum of 
WISE J084502.48-545808.4, candidate counterpart of 
2FGL J0844.8-5459, show an absorption feature {Xobs 
= 6364 A with EWobs = 2.2 A) that we are not able to 
clearly identify (see Figure 15). We classify the WISE 
J162638.15-763855.5 counterpart to 2FGL J1626.0- 


7636 as a BL Lac since its emission lines have EW < 5 
A. The detection of emission ([O I] with EWobs = 2.3 
A, the [S II] doublet AAobs = 7421 - 7438 A with 
EWobs = 2.5 - 2.4 A) and absorption lines (G band. 
Mg I with EWobs = 4.9 A and Na I with EWobs = 3.3 
A), enable us to measure a redshift of 0.1050 (see 
Figure 20). The optical spectra of the AGUs associ¬ 
ated with the y-ray source 2FGL J1849.5-r2744 are 
also published in Shaw et al. (2013a). The source 
is listed as a BZB in the 2LAC catalog, although 
the optical spectrum was not yet available. The op¬ 
tical spectrum collected by us for this object (WISE 
J18493L74-r274800.8) is nearly featureless (see Fig¬ 
ure 21), so we are not able to confirm the lower limit z 
estimate (i.e., 1.466) of this source proposed by Shaw 
et al. (2013a) on the basis of Mg II absorption dou¬ 
blet. Thanks to our optical spectrum of the WISE 
J195945.66-472519.3, candidate counterpart to 2FGL 
J1959.9-4727, we can put a lower limit on its redshift 
of 0.519 based on the detection of an intervening dou¬ 
blet system of Mg II (AAobs = 4246 - 4256 A with 
EWobs = 1.4 - 0.9 A; see also Figure 22). 

4.3. BL Lac objects 

Details for the BL Lacs in our sample are listed be¬ 
low. Table 1 reports their Roma-BZCAT name and 
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the name of the WISE counterpart with the coordi¬ 
nates. Multifrequency notes are not presented in this 
table since they are already extensively discussed in 
the Roma-BZCAT. All the BZBs that are also de¬ 
tected by Fermi belong to the sample named: locus 
(i.e. the three-dimensional region in the WISE color 
space occupied by the y-ray-emitting blazars, well 
separated from the regions occupied by other classes, 
D’Abrusco et al. 2013), with the only exception being 
BZB J2039-I-5219. They were used in D’Abrusco et 
al. (2013) to build the method to search for blazar- 
like sources within the positional uncertainty region 
of the Fermi UGSs. Thus all their WISE counterparts 
have the IR colors consistent with the Fermi blazar 
population. There are five objects listed in Table 1 
as BL Lac candidates for which no optical spectra are 
present in literature that allowed a firm classification. 
Our spectroscopic observations confirmed that all the 
BL Lac candidates have featureless optical spectra as 
shown in Ligures 23-30. In addition, the remaining 
three objects, namely BZB J1217-r3007, BZB J1359- 
3746 and BZB J1553-3118, that were indeed classi¬ 
fied BZBs with an uncertain redshift estimate when 
the Roma-BZCAT was released. Lor one BZB, BZB 
J0244-5819, we have been able to estimate a lower 
limit on its redshift, z>0.265, on the basis of absorp¬ 
tion features in its optical spectrum. Our observation 
of WISE J024440.30-581954.5 show the doublet Ca 
H-rK {AAobs = 4975 - 5017 A with EWobs = 0.7-1 A), 
G band. Mg I {EWobs = 1.5 A) and Na I {EWobs = 1.1 
A) absorption features that could be due to the host 
galaxy and/or to intervening systems (see Ligure 23). 
We note that BZB J1359-3746, BZB J1649-r5235 and 
BZB J1809-r2910 also have published optical spectra 
for their candidate counterparts released in Shaw et al. 
(2013b). In addition, Shaw et al. (2013b) reported 
the spectrum of BZB J1359-3746, with the detection 
of the Calcium break that allowed them to determine a 
redshift of 0.334. We cannot confirm this result since 
our optical spectrum of WISE J135949.71-374600.7 
shows a featureless continuum (see Ligure 25). We 
also present the spectrum of the low-energy counter¬ 
part of BZB J1553-3118, which was also published 
in Masetti et al. (2013) with a redshift of 0.132. In 
our lower signal to noise observation, its IR coun¬ 
terpart WISE J135949.71-374600.7 appears feature¬ 
less, as shown in Ligure 26. The optical spectrum 
of WISE J203923.51-r521950.1, candidate counterpart 
of BZB J2039-r5219, show an absorption feature at 
Aobs = 6212 A with EWobs = 4.8 A. We do not have 


a clear identification for this feature (see also Ligure 
30). 

5. Summary and conclusions 

We present here the results of our 2014 spectro¬ 
scopic campaign carried out in the Northern hemi¬ 
sphere with the KPNO and Palomar telescopes, and in 
the Southern hemisphere with the SOAR and Magellan 
telescopes. The main goal of our campaign is to con¬ 
firm the nature of sources selected for having IR colors 
or low radio frequency spectra (i.e., below ~lGHz) 
similar to known Eermi-doitcitd blazars and lying 
within the positional uncertainty regions of the UGSs 
through optical spectroscopic observations. Given the 
positional uncertainty of the UGSs, there could be 
a possible contamination by different source classes 
in these selection procedures of y-ray blazar can¬ 
didates (e.g.. Stern & Assef 2013) and spectroscopic 
observations are the only way to remove such de¬ 
generacy (e.g., Masetti et al. 2013; Shaw et al. 2013a; 

Shaw et al. 2013b; Paggi et al. 2014; Massaro et al. 2014a) 
and discover new gamma-ray blazars. The confir¬ 
mation of the blazar nature of these selected ob¬ 
jects will improve/refine future associations for the 
Fermi catalog and will also yield to understand 
the efficiency and completeness of the association 
method based on the WISE colors once our cam¬ 
paign is completed. Our spectroscopic observations 
could potentially allow us to obtain redshift esti¬ 
mates for the UGS candidate counterparts. During 
our campaign we also observed several AGUs, as 
defined according to the Fermi catalogs (see e.g., 
Ackermann et al. 2011a; Nolan et al. 2012), to verify 
if they are indeed blazars. In addition we observed sev¬ 
eral sources that already belong to the Roma-BZCAT 
but were classified as BL Lac candidates due to the 
lack of optical spectra available in literature, or were 
BL Lac objects with uncertain redshift estimates, when 
the Roma-BZCAT v4.1 was released. 

We observed a total of 27 targets. The results of this 
spectroscopic campaign are reported as follows: 

• In the sample of potential counterparts for 
the UGS, selected based on their IR colors 
(Massaro et al. 2011a; D’Abrusco et al. 2013; 
Massaro et al. 2013a) and on the basis of their 
fiat radio spectra below ~ 1 GHz (Massaro et al. 2013b; 
Nori et al. 2014; Massaro et al. 2014b), we con¬ 
firm the blazar nature of 8 of 9 UGSs. Among 
them, six are clearly BL Lacs presenting fea- 
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tureless optical spectra. The remaining two 
are QSOs. One potential counterparts for 
2FGL J1745.6+0203 found on the basis of 
the IR color selection method, namely WISE 
J174526.95+020532.6, is a QSO that probably 
contaminates our selection method. 

• We classify WISE J173052.85-035247.2 as a BL 
Lac although this source is not associated with 
any NVSS counterpart, as expected for BL Lac 
objects. The detection of an intervening doublet 
system of Mg II enabled us to set a lower limit 
on its redshift of 0.776. 

• All the sources that belong to our AGU sub¬ 
sample have a blazar nature. Two of them are 
QSOs, while for one of them, namely WISE 
J162638.15-763855.5, we have been able to de¬ 
tect some features in the optical spectrum, both 
in emission and in absorption, leading to a red- 
shift measurement of 0.1050. The AGU asso¬ 
ciated with the W/5'£’ J18493L74+274800.8 has 
been classified as a BZB in the 2LAC paper but 
at that time there was no optical spectra. Our ob¬ 
servations confirm its BL Lac nature. We have 
been also able to set a lower limit for the AGU 
associated with WISE J195945.66-472519.3 of 
0.519 thanks to the detection in its optical spec¬ 
trum of a Mg II intervening system. 

• Within the Roma-BZCAT sources we found five 
BL Lac candidates that thanks to the collected 
optical spectra are all confirmed BL Lacs. For 
one of them, BZB J0244-5819, we have been 
able to set a lower limit on their redshifts on the 
basis of absorption systems (Ca H+K, G band. 
Mg I and Na I) that could be due to the host 
galaxy and/or to intervening systems. For the re¬ 
maining four BZBs listed in the Roma-BZCAT 
with uncertain redshift estimates, we were not 
able to obtain any redshift values. 
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Fig. 1.— Optical images, 5' on a side, of 9 of the WISE sources selected in this paper for optical spectroscopic follow¬ 
up. The proposed optical counterparts are indicated with red marks and the fields are extracted from the DSS-II-Red 
survey. Object name, image scale and orientation are also reported in each panel. 
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Fig. 2.— Same as Figure 1 . 
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Fig. 3.— Same as Figure 1. 
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Fig. 4.— Upper panel: optical spectra observed at 
SOAR of WISE J112325.37-252857.0, potential coun¬ 
terpart associated with 2FGL J1123.3-2527, classified 
as a QSO at z=0.148 on the basis of the identification 
of the emission lines shown in the spectra. The aver¬ 
age SNR is also reported in the figure. Lower panel: 
the normalised spectrum is shown. The symbol 0 in¬ 
dicates atmospheric telluric features. 
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Fig. 5.— Upper panel: optical spectra observed at 
SOAR of W/5'F’ J125949.80-374858.1, potential coun¬ 
terpart associated with FGL J1259.8-3749, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 6.— Upper panel: optical spectra observed at 
SOAR of J132840.61-472749.2, potential coun¬ 
terpart associated with 2FGL J1328.5-4728, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 8.— Upper panel: optical spectra observed at 
SOAR of J134706.88-295842.4, potential coun¬ 
terpart associated with 2FGL J1347.0-2956, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 7.— Upper panel: optical spectra observed at 
KPNO of J134042.01-041006.8, potential coun¬ 
terpart associated with 2FGL J1340.5-0412, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 


Fig. 9.— Upper panel: optical spectra observed at 
KPNO of J173052.85-035247.2, potential coun¬ 
terpart associated with 2FGL J1730.6-0353, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 10.— Upper panel: optical spectra observed 
at KPNO of WISE J174507.82+015442.4, poten¬ 
tial counterpart associated with 2FGL J1745.6-r0203, 
classified as a QSO at z=0.078 on the basis of the iden¬ 
tification of the emission lines visible in the spectra. 
Lower panel: as in Figure 4. 
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Fig. 11.— Upper panel: optical spectra observed 
at KPNO of WISE J174526.95-r020532.6, poten¬ 
tial counterpart associated with 2FGL J1745.6-h0203, 
classified as a QSO at z=0.335 on the basis of the iden¬ 
tification of the emission lines visible in the spectra. 
Lower panel: as in Figure 4. 
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Fig. 12.— Upper panel: optical spectra observed 
at KPNO of WISE J183535.34-hl34848.8, poten¬ 
tial counterpart associated with 2FGL J1835.4-h 1349, 
classified as a BL Lac on the basis of its featureless 
continuum. Lower panel: as in Figure 4. 
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Fig. 13.— Upper panel: optical spectra observed 
at Palomar of WISE J025333.64-h321720.8, poten¬ 
tial counterpart associated with 2FGL J0253.4-h3218, 
classified as a QSO at z=0.859 on the basis of the iden¬ 
tification of the emission lines listed on the spectra. 
Red and blue spectra show the two sides of the dual¬ 
beam spectrograph. Lower panel: as in Figure 4. 
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Fig. 14.— Upper panel: optical spectra observed at 
SOAR of WISE J074642.31-475455.0, potential coun¬ 
terpart associated with 2FGL J0746.5-4758, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 



Fig. 15.— Upper panel: optical spectra observed at 
SOAR of WISE J084502.48-545808.4, potential coun¬ 
terpart associated with 2FGL J0844.8-5459, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 16.— Upper panel: optical spectra observed 
at Palomar of WISE J085654.85-r714623.8, poten¬ 
tial counterpart associated with 2FGL J0856.0-r7136, 
classified as a QSO at z=0.542 on the basis of the iden¬ 
tification of the emission lines visible in the spectra. 
Red and blue spectra show the two sides of the dual¬ 
beam spectrograph. Lower panel: as in Figure 4. 



Fig. 17.— Upper panel: optical spectra observed at 
SOAR of WS^Jl 14141.80-140754.6, potential coun¬ 
terpart associated with 2FGL J1141.7-1404, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 18.— Upper panel: optical spectra observed at 
SOAR of lF/5'£’ J123824.39-195913.8, potential coun¬ 
terpart associated with 2FGL J1238.1-1953 classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 20.— Upper panel: optical spectra observed at 
SOAR of J162638.15-763855.5, potential coun¬ 
terpart associated with 2FGL J1626.0-7636, classified 
as a BL Lac since its emission and absorption lines 
have EW ] 5 A. The detection of emission ([O I] and 
the [S II] doublet) and absorption lines (G band, Mg 
I and Na I) has enabled us to estimate a redshift of 
0.1050. Lower panel: as in Figure 4. 
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Fig. 19.— Upper panel: optical spectra observed at 
SOAR of WS^J140609.60-250809.2, potential coun¬ 
terpart associated with 2FGL J1406.2-2510, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 


Fig. 21.— Upper panel: optical spectra observed 
at KPNO of WISE J18493L74-r274800.8, poten¬ 
tial counterpart associated with 2FGL J1849.5-r2744, 
classified as a BL Lac on the basis of its featureless 
continuum. Lower panel: as in Figure 4. 
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Fig. 22.— Upper panel: optical spectra observed at 
SOAR of J195945.66-472519.3, potential coun¬ 
terpart associated with 2FGL J1959.9-4727, classified 
as a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 23.— Upper panel: optical spectra observed 
at Magellan of WISE J024440.30-581954.5, potential 
counterpart associated with BZB J0244-5819, classi¬ 
fied as a BL Lac since its absorption lines have EW j 5 
A. The detection of absorption lines (Ca H-hK, G band, 
Mg I and Na I) has enabled us to estimate a lower limit 
on its redshift of 0.265. Lower panel: as in Figure 4. 



Fig. 24.— Upper panel: optical spectra observed 
at KPNO of WISE J121752.08-r300700.6, potential 
counterpart associated with BZB J1217-r3007, classi¬ 
fied as a BL Lac on the basis of its featureless contin¬ 
uum. Lower panel: as in Figure 4. 
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Fig. 25.— Upper panel: optical spectra observed at 
SOAR of J135949.71-374600.7, potential coun¬ 
terpart associated with BZB J1359-3746, classified as 
a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 
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Fig. 26.— Upper panel: optical spectra observed at 
SOAR of lF/5'£’ J155333.56-311830.9, potential coun¬ 
terpart associated with BZB J1553-3118, classified as 
a BL Lac on the basis of its featureless continuum. 
Lower panel: as in Figure 4. 



wavelength [A] 

Fig. 28.— Upper panel: optical spectra observed 
at KPNO of WISE J170209.63-r264314.7, potential 
counterpart associated with BZB J1702-h2643, classi¬ 
fied as a BL Lac on the basis of its featureless contin¬ 
uum. Lower panel: as in Figure 4. 
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Fig. 27.— Upper panel: optical spectra observed 
at KPNO of WISE J164924.98-r523515.0, potential 
counterpart associated with BZB J1649+5235, classi¬ 
fied as a BL Lac on the basis of its featureless contin¬ 
uum. Lower panel: as in Figure 4. 



Fig. 29.— Upper panel: optical spectra observed 
at KPNO of WISE J180945.39+291019.8, potential 
counterpart associated with BZB J1809+2910, classi¬ 
fied as a BL Lac on the basis of its featureless contin¬ 
uum. Lower panel: as in Figure 4. 
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Fig. 30.— Upper panel: optical spectra observed 
at KPNO of WISE J203923.51+521950.1, potential 
counterpart associated with BZB J2039+5219, classi¬ 
fied as a BL Lac on the basis of its featureless contin¬ 
uum. Lower panel: as in Figure 4. 
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